FLUCTUATION ANALYSIS FOR THE LOSS FROM DEFAULT 
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Abstract. We analyze the fluctuation of the loss from default around its large portfolio limit in a class of 
reduced-form models of correlated firm-by-firm default timing. We prove a weak convergence result for the 
fluctuation process and use it for developing a conditionally Gaussian approximation to the loss distribution. 
Numerical results illustrate the accuracy and computational efficiency of the approximation. 



1. Introduction 

Reduced-form point process models of correlated firm-by-firm default timing arc widely used to measure 
the credit risk in portfolios of defaultable assets such as loans and corporate bonds. In these models, 
defaults arrive at intensities governed by a given system of stochastic differential equations. Computing the 
distribution of the loss from default in these models is often challenging. Transform methods are tractable 
only in special cases, for example, if defaults are assumed to be conditionally independent. Monte Carlo 
simulation methods have a much wider scope but can be slow for the large portfolios and longer time 
horizons common in industry practice. A major US bank might easily have 20, 000 wholesale loans, 50, 000 
to 100, 000 mid-market and commercial loans, and derivatives trades with 10, 000 to 20, 000 different legal 
entities. Simulation of such large pools is extremely burdensome. 

This paper develops a conditionally Gaussian approximation to the distribution of the loss from default 
in large pools. The approximation applies to a broad class of empirically motivated reduced- form models in 
which a name defaults at an intensity that is influenced by an idiosyncratic risk factor process, a systematic 
risk factor process X common to all names in the pool, and the portfolio loss. It is based on an analysis 
of the fluctuation of the loss around its large portfolio limit, i.e., a central limit theorem (CLT). More 
precisely, we show that the signed-measurc-valued process describing the fluctuation of the loss around its 
law of large numbers (LLN) limit weakly converges to a unique distribution-valued process in a certain 
weighted Sobolcv space. The limiting fluctuation process satisfies a stochastic evolution equation driven by 
the Brownian motion governing the systematic risk X and a distribution-valued martingale that is centered 
Gaussian given X. The fluctuation limit, and thus the resulting approximation to the loss distribution, is 
Gaussian only in the special case that the names in the pool are not sensitive to X . In the general case, the 
approximation is conditionally Gaussian. 

The weak convergence result proven in this paper extends the LLNs for the portfolio loss established 
in our earlier work Giesecke, Spiliopoulos & Sowers (2013) and Giesecke, Spiliopoulos, Sowers & Sirignano 
(2012). The fluctuation analysis performed here involves challenging topological issues that do not arise in 
the analysis of the LLN. Firstly, the fluctuations process takes values in the space of signed measures. This 
space is not well suited for the study of convergence in distribution; the weak topology, which is the natural 
topology to consider here, is in general not mctrizable (see, for example, Del Barrio, Deheuvels & Geer 
(2007)). We address this issue by analyzing the convergence of the fluctuations as a process taking values in 
the dual of a suitable Sobolev space of test functions. Weights are introduced in order to control the norm 
of the fluctuations and establish tightness. Related ideas appear in Metivier (1985), Fernandez & Melcard 
(1997), Kurtz & Xiong (2004) and other articles, but for other systems and under different assumptions. 
Additional issues are the growth and degeneracies of the coefficients of our stochastic system, which make 
it difficult to identify the appropriate weights for the Sobolev space in which the convergence is established 
and uniqueness of the limiting fluctuation process is proved. Purtukhia (1984), Gyongy & Krylov (1990), 
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Krylov & Lototsky (1998), Lototsky (2001), Kim (2009), and others face similar issues in settings that are 
different from ours. The approach we pursue is inspired by the class of weights that were introduced by 
Purtukhia (1984) and further analyzed by Gydngy & Krylov (1990). 

The fluctuation analysis leads to a second-order approximation to the loss distribution that can be sig- 
nificantly more accurate than an approximation obtained from just the LLN. Our numerical results, which 
are based on a method of moments for solving the stochastic evolution equation governing the fluctuation 
limit, confirm that the second-order approximation is much more accurate. The second-order approximation 
is even accurate for relatively small portfolios in the order of hundreds of names or when the influence of the 
systematic risk process X is relatively low. The second-order approximation also improves the accuracy in 
the tails. The effort of computing the second-order approximation exceeds that of computing the first-order 
approximation, but is still much lower than compared to directly simulating the pool. 

Prior research has established various weak convergence results for interacting particle systems represented 
by reduced-form models of correlated default timing. Dai Pra, Runggaldier, Sartori & Tolotti (2009) and 
Dai Pra & Tolotti (2009) study mean-field models in which an intensity is a function of the portfolio loss. In a 
model with local interaction, Giesecke & Weber (2006) take the intensity of a name as a function of the state 
of the names in a specified neighborhood of that name. In these formulations, the impact of a default on the 
dynamics of the surviving constituent names, a contagion effect, is permanent. In our mean-field model, an 
intensity depends on the path of the portfolio loss. Therefore, the impact of a default may be transient and 
fade away with time. There is a recovery effect. Moreover, the system which we analyze includes firm-specific 
sources of default risk and addresses an additional source of default clustering, namely the exposure of a 
name to a systematic risk factor process. This exposure generates a random limiting behavior for the LLN 
and leads to a fluctuation limit which is only conditionally Gaussian. 

There are several other related articles. Cvitanic, Ma & Zhang (2012) prove a LLN for a mean-field system 
with permanent default impact, taking an intensity as a function of an idiosyncratic risk factor, a systematic 
risk factor, and the portfolio loss rate. The risk factors follow diffusion processes whose coefficients may 
depend on the portfolio loss. Bush, Hambly, Haworth, Jin & Reisinger (2011) establish a LLN for a system 
in which defaults occur at first hitting times of correlated diffusion processes. Conditional on a correlating 
systematic risk factor governed by a Brownian motion, defaults occur independently of one another. Gamier, 
Papanicolaou & Yang (2012) analyze large deviations in a system of diffusion processes that interact through 
their empirical mean and have a stabilizing force acting on each of them. Fouque & Ichiba (2012) study 
defaults and stability in an interacting diffusion model of inter-bank lending. 

The rest of this paper is organized as follows. Section [5] describes the class of reduced-form models of 
correlated default timing which we analyze. Section [3] reviews the law of large numbers proved by Giesecke 
et al. (2012) for the portfolio loss in these models. Section [4] states our main weak convergence result for the 
fluctuation process, essentially a central limit theorem for the loss. Section [5] provides a numerical method 
for solving the stochastic evolution equation governing the fluctuation limit and provides numerical results. 
Sections [BE] are devoted to the proof of the main result. There is an appendix. 



2. Model and Assumptions 

We analyze a class of reduced-form point process models of correlated default timing in a pool of firms 
("names"). Models of this type have been studied by Cvitanic et al. (2012), Dai Pra et al. (2009), Dai Pra 
& Tolotti (2009), Giesecke et al. (2013), Giesecke et al. (2012), and others. In these models, a default is 
governed by the first jump of a point process. The point process intensities follow a system of SDEs. 

Fix a probability space (fi,^,P). Let {VF n }„<=N be a countable collection of independent standard 
Brownian motions. Let {e„} nS N be an i.i.d. collection of standard exponential random variables which are 
independent of the W n 's. Finally, let V be a standard Brownian motion which is independent of the W n, s and 
e„'s. Each W n will represent a source of risk which is idiosyncratic to a specific name. Each e n will represent 
a normalized default time for a specific name. The process V will drive a systematic risk factor process 
to which all names are exposed. Define V t = cr (V s , < s < t) V Af and J? n = a {(V a ,W?), 0<s<t)VAf, 
where Af contains the P-null sets. Lastly, we will denote by FY the conditional law given Y, where Y may 
represent Vt, for example. 



Fix N € N, n € {1, 2, . . . , N} and consider the following system: 
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Here, x ' 1S the indicator function. The initial value x of X is hxed. The ajy,n, Ajv.m <7iV,nj /3^ „, /?jv „, A 0l Ar,n 
are parameters that are specified below. The process L N represents the loss rate in a portfolio of N names, 
assuming a loss at default of one unit. The process A^' n represents the intensity, or conditional default rate, 
of the n-th name in the pool. More precisely, A"'" is the density of the Doob-Mcyer compensator to the 
default indicator 1 — M t '", where M t '" = X{r N > n >t}- That is, a martingale is given by M t ' n + / X^ ,n M^' n ds. 
The results in Section 3 of Giesecke et al. (2013) imply that the system (JTJ) has a unique solution such that 
Af '" > for every N € N, n G {1, 2, . . . , N} and t > 0. Thus, the model is well-posed. 

The default timing model ([1]) addresses several channels of default clustering. An intensity is influenced 
by an idiosyncratic source of risk represented by a Brownian motion W n , and a source of systematic risk 
common to all firms-the diffusion process X. Movements in X cause correlated changes in firms' intensities 
and thus provide a channel for default clustering emphasized by Das, Duffie, Kapadia & Saita (2007) for 
corporate defaults in the U.S. The sensitivity of X N ' n to changes in X is measured by the parameter /3jy n e K. 
The second channel for default clustering is modeled through the feedback ("contagion") term ^ n dLf . 
A default causes a jump of size fl% n /N in the intensity A^'™, where fi% n € M + = [0,oo). Due to the 
mean- reversion of A '", the impact of a default fades away with time, exponentially with rate aN,n € K+- 
Azizpour, Giesecke & Schwenkler (2010) have found self-exciting effects of this type to be an important 
channel for the clustering of defaults in the U.S., over and above any clustering caused by the exposure of 
firms to systematic risk factors. Giesecke & Schwenkler (2011) develop and analyze likelihood estimators of 
the parameters of point process models such as ([1]). 

We allow for a heterogeneous pool; the intensity dynamics of each name can be different. We capture 
these different dynamics by defining the parameter "types" 

(2) P '" = (aN,n.,XN,n,VN,n,PN,nipN,n)- 

The p-^^'s take values in parameter space V = Ri x R. For each iVeN, define 

~N,n def ,N,n \N,n\ 

Pt = (P >A t ) 

for all n € {1,2,..., N} and t > 0. Define P = Px R + . The vector p = (pj, . . . , p^) represents a random 
environment that addresses the heterogeneity of the system. 

Condition 2.1. We assume that the Pg are i.i.d. random variables with common law v and that v has 
compact support in V . In particular, and to avoid inessential complications, we assume that all elements of 
the random vector pg are bounded from above and from below (for those that take values in RJ by some K , 
for all n 6 {1, 2, . . . , N}. Moreover, we assume that {pg} is independent of {W n }, {c n } and V . 

This formulation of heterogeneity generalizes that of Giesecke ct al. (2012). They assume that n = 
limjv^oo jj E n =i ^p"-" and A o = limAT^oo i J2 n =i ^ 0]N ,„ exist (in £?(P) and ^(R+), respectively). Their 
formulation is a special case of the one proposed here, where the law v = tt X A . 

Condition 2.2. We assume that Ej*[\b (X s )\ 2 + |cr (A" s )| 4 ] ds < oo for allt>0. 

3. Law of large numbers 

Except for special cases of little practical interest, the distribution of the portfolio loss L N in the system 
([T]) is difficult to compute. We are interested in an approximation to this distribution for the large portfolios 
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common in practice, i.e., for the case that TV is large. Giesecke et al. (2012) prove a law of large numbers 
(LLN) for the loss in the system ([1]) and use it for developing a first-order approximation. We first review 
this result and then extend it in Section U by analyzing the fluctuations of the loss around its LLN limit. 
The fluctuation analysis will allow us to construct a more accurate second-order approximation. 
To outline the LLN, define 

i N 

JV def J- v— v r , JV.n 

Lit = > 0~N,nM+ ' ; 

71=1 

this is the empirical distribution of the type and intensity for those names which are still "alive." We note 
that n^ is a sub-probability measure. Since 

(3) Lf = i- M f(n 

it suffices to study the limiting behavior of the measure-valued process {/xf ,t G [0,T]}jvgn for some fixed 
horizon T > 0. Let E be the collection of sub-probability measures (i.e., defective probability measures) on 
V) i.e., E consists of those Borel measures v on V such that v(V) < 1. Topologizing E in the usual way (by 
projecting onto the one-point compactification of V\ see Chapter 9.5 of Royden (1988)) we obtain that E a 
Polish space. Thus, fi N is an element of De[0, oo) where D is the Skorokhod space (i.e., De[0, oo) is the set 
of RCLL processes on [0, oo) taking values in E). 

Further, for p = (p, A) where p = (a, \,a, f3 c , (3 s ) G V and / G C 2 (V) (the space of twice continuously 
diffcrentiable, bounded functions), define the operators 

(4) (A/)( P ) = 2 ct2a 9a4(p) - "( A - ^(p) - a /(p) 

(5) 0C2/)(p) = |£(P) 

(6) mm) = (3 s \b (x)^(p) + -(/3 s ) 2 A 2 a 2 (x)^(p) 

(7) (£lf)(p)=p s \a a (x)^.(p). 

Also define 

Qi(p) d =A/3 c . 

The generator C\ corresponds to the diffusive part of the intensity with killing rate A, and £2 is the macro- 
scopic effect of contagion on the surviving intensities at any given time. The operators £| and L\ are related 
to the exogenous systematic risk X. For a measure v G Db[0,oo), we also specify the inner product 



</,!/>= / /(p)<Mp). 

Jv 

The law of large numbers of Giesecke ct al. (2012) states that /uf weakly converges to /2 t in De[0,T]. To 
rigorously formulate the result we need to use the weak form (see also Lemma I8.3|) . In particular, for all 
/ G C'friV), the evolution of ft. is governed by the measure evolution equation 

(8) d(f,fi t ) E ^[(C 1 f,fit) E + (Qi,iit) E (C2f,iit) E + (cf t f,fk) E }dt + ^f,fk) E dV t , a.s. 

The LLN suggests an approximation to the distribution of the loss L in large pools by the large pool limit: 

(9) L»~L t = l-ji t {f). 



4. Main Result: Fluctuations Theorem 

In order to improve the first-order approximation Q, we analyze the fluctuations of ll n around its large 
pool limit ft. As is indicated by the proof of Lemmas 18.41 and 18.51 (see also Giesecke et al. (2012)), for an 
appropriate metric, the sequence {v/V (fi^ — jm) '■ TV < 00} is stochastically bounded. Hence, it is reasonable 
to define the scaled fluctuation process E N by 

(10) Sf = y/N(tf - fit). 



In Theorem 14.11 below we will show that the signed-measure-valued process S^ weakly converges to a 
fluctuation limit process S in an appropriate space. 

The analysis of the limiting behavior of the fluctuation process ([TT)|) involves issues that do not occur 
in the treatment of the LLN. In particular, even though the fluctuation process is a signed-measure-valued 
process, its limit process {S.} is distribution- valued in an appropriate space. The space of signed measures 
endowed with the weak topology is not metrizablcu The difficulty is then to identify a rich enough space, 
where tightness and uniqueness can be proven. It turns out that we have to consider the convergence in 
weighted Sobolev spaces. Here, several technical challenges arise. These are mainly due to the growth and 
degeneracies of the coefficients of the system (TTJ), which make it difficult to identify the correct weights. The 
spaces that we consider are Hilbert spaces. For the sake of clarity of presentation the appropriate Hilbert 
spaces will be defined in detail in Section [7] For the moment, we mention that the space in question is 
denoted by Wq(w, p), with w and p the appropriate weight functions, JeN and W -J (iu, p) will be its dual. 
A precise definition is given in Section [Jj 

We need to introduce several additional operators to state our weak convergence result. For p = (p, A) e V 
and / e C%(T), we define 

(ii) (& iM /)(p) = (A/)(p) + mm) + (Qi,n) E (£ 2 /)(p) + <£ 2 /, v) E Ci(p) 

(12) 

(13) (£ 6 (/, 3 ))(p) = /(p).g(p)A 

(14) (£ 7 /)(p) = /3 c /(p)A 

(15) fi 2 (p) = A(/3 C ) 2 . 

The main result of this paper is the following theorem. 

Theorem 4.1. Let J = 2. The sequence {S^,i G [0,T]}jv e N is relatively compact in D w -(j+2), J0,T]. 
For any subsequence of this sequence, there exists a subsubsequence that converges in distribution with limit 
{S t ,i G [0,T]}. Any accumulation point S satisfies the stochastic evolution equation 

(16) </,St) = </,Ho>+ / (Gx s jiJ,Z s )ds + J (c?'f,E a )dV a + (f,M t ), a.s. 

for any f G W (w,p), where M. is a distribution-valued martingale with predictable variation process 

(17) \{f,M)]t= I \(C 5 (fJ),p s ) + (£ 6 (/,/),/2 s ) + (C 2 f,p s ) 2 (Q 2 ,p s ) -2(£ 7 /,/2.) (£a/,jS.) 



r 2 d f f^ d 9 , 



(£ 5 (/, 3 ))(p) = (T 2 ^-(p)^(p)A 



ds. 

Moreover, conditional on the a-algebra Vt, M-t is centered Gaussian with covariance function, for f,g G 



Wq (w,p), given by 



Cov 



(18) 



(f,M tl ),(g,M t2 ) 



V, 



*lVt 2 



= E 



t 1 M 2 



[(£5{f,g),Vs) + (£6(f,g),Hs) + (C 2 f,n a ) (C 2 g,ti s ) (0.2,1*8 



(£ 7 g, p s ) (£ 2 /, p s ) - (C 7 f, p s ) (£23, Ps)} ds V tl vt 2 



(J+2)/ 



Finally, the limiting stochastic evolution equation A 16]) has a unique solution in W (w,p) and thus the 

limit accumulation point E. is unique. 



In general, a topological space is Polish if and only if its topology can be defined by a metric for which it is complete 
and separable. A locally compact space is metrizable if and only if it has a countable base in which it is Polish. The weak 
convergence of probability measures (i.e., finite non-negative measures) on a separable metric space can be defined by the 
Prokhorov metric. Thus, the weak topology on the space of probability measures on a locally compact space E is metrizable 
if and only if E is Polish. However, issues arise if one replaces the space of finite non-negative measures with the space of 
finite signed measures. The space of bounded and continuous functions on E, endowed with the sup-norm is a Banach space. 
Its topological dual endowed with the weak* -topology coincides with the set of finite signed measures on E endowed with 
the weak topology. However, the topological dual of an infinite dimensional Banach space is not metrizable, even though, by 
the Banach-Alaoglu theorem, any weak*-compact subset of such a topological dual will be metrizable. For a more thorough 
discussion of these issues, sec Del Barrio et al. (2007), Remark 1.2. 
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Remark 4.2. Clearly, when P% n = for all N,n, then the limiting distribution-valued martingale M. is 
centered Gaussian with covariance operator given by the (now deterministic) term within the expectation in 
H8\). Also, we remark that the operators G x ,ft an d C% in the stochastic evolution equation 1116]) are linear. 
Conditionally on the systematic risk X , equation 111 6]) is linear. 

The proof of Theorem 14.11 is developed in Sections [5] through [HI In Section [SI we identify the limiting 
equation and prove the convergence theorem based on the tightness and uniqueness results of Sections [8] and 
Section [9) In Section [3 we discuss the Sobolev spaces we are using and the necessity of introducing weights. 
In Section|8l we prove that the family {S t ' n ,t £ [0,T]}jv 6 n is relatively compact in D _(j +2 ), . [0,T], see 

Lemma \8M In Section [9l we prove uniqueness of (fl~6)) in W (w,p), see Theorem 19.71 

The fluctuation analysis leads to a second-order approximation to the distribution of the portfolio loss 
L N in large pools. The weak convergence established in Theorem 14.11 implies that 

(19) F(VN(L? -L t )>i)^P(EtCP)<-£) 

for large N. Theorem 14. II and its proof in Section [5] imply that (H , V, p,) weakly converges to (S, V, p,). This 
motivates the approximation 

(20) "f = ^ +ft ^ +fc 



which implies a second-order approximation for the portfolio loss: 
(21) L?&L t --j=E t (P). 



The next section develops, implements and tests a numerical method for computing the distribution of 
L t — E. t (V)/vN, and numerically demonstrates the accuracy of the approximation (|21j) . 

5. Numerical Solution 

The numerical solution of the stochastic evolution equation for the fluctuation limit is not standard. There 
surprisingly exists little literature on solving this class of problems. Solutions do not exist in L 2 , so standard 
Galerkin methods are not applicableLi 

5.1. Method of Moments. We provide a method of moments for solving for the fluctuation limit S('P). 
Along with the LLN limit p,(P) = 1 — L, the fluctuation limit yields the approximation ([21]) . The method 
extends the numerical approach for computing p(V) developed by Giesecke et al. (2012). We first present it 
for a homogeneous portfolio; i.e., set v = <5p where po 6?. Write S t (<iA) for the solution of the stochastic 
evolution equation (|16[) ; for notational convenience, we do not explicitly show the dependence on the fixed 
parameters. For k € N, we define the "fluctuation moments" 

/•OO 

(22) v k {t) = / X k E t (d\). 

Jo 

We are interested in vo(i), which is equal to the fluctuation limit H((R+). We also define the moments Uk(t) 
of the solution p t (d\) of the stochastic evolution equation ([8]). again not explicitly showing the dependence 
on the fixed parameters: 

/•OO 

(23) u k (t) = / \ k jLt{d\). 



An example of the unique challenges posed by these equations is the simple case of i.i.d. Brownian motions starting from zero 
on the real line. The fluctuation limit Ct{ x ) f° r this system satisfies the stochastic evolution equation rf(/, £t) = ™ (/"i C')^* + 
d(f, Wt), where Wt is a Gaussian process with covariance Cov[(/, Wt), (g, W s )] = f^ At J R f'(x)g'(x)(2Tru)- 1 / 2 e- x2 / < - 2u Uxdu. 

Suppose there exists a solution £(£, x) which satisfies the SPDE df = ~ fdt-\-d\Vf Then, C,t(x) = C(*! x)dx and one could solve 
the stochastic evolution equation by solving this SPDE. However, challenges are immediately evident. The Green's function 
solution Q(t,x) = f J m dt'dx'G(t' ,x';t,x)dW t i to the SPDE has infinite variance, indicating that the Green's solution is not 
in fact a solution to this SPDE. It is interesting to note that the lack of a Green's function solution stems directly from the 
covariance structure; if one had f(x)g(x) instead of f'(x)g'(x), there would be a finite variance Green's function solution. 
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The zero-th moment gives the limiting loss L since L t = 1 — Uo(t). Both the zero-th LLN moment uo(t) 
and the zero-th fluctuation moment i>o(t) arc required to compute the approximation (|2 1 [) . The moments 
{ u k(t)}kLo satisfy a system of SDEs, 

du k (t) = {«*(*)(- ak + l3 s b Q {X t )k + 0.5((3 s ) 2 a 2 {X t )k(k - 1)) 

(24) + u k -i{t){Q.ba 2 k{k - 1) + aA/c + /3 c fc Wl (t)) - u k+1 {t)}dt + P s a (X t )ku k {t)dV t , 

Ufe(O) = / X k p {dX), 

Jo 

see Giesecke et al. (2012). The fluctuation moments {v k (t)} k x L 1 can also be shown to satisfy a system of 
SDEs using the stochastic evolution equation (|16p . Taking a test function f = X k in (fT6|) . we have 

cfe fc (£) = /3 C ku k -i(t)vi(t)dt + [0.5cr 2 /c(/c - 1) + aXk + k/3 C Ui(t)] v k -i(t)dt - v k+l {t)dt 

(25) + [kp s b (X t ) -ka + 0.5k(k - 1) (p s a (X t )) 2 } v k (t)dt + k/3 s a (X t )v k {t)dV t + dM k {t), 

v k (0) = / X k E (dX). 



J 
where M. k (t) = (X k ,M.t) an d its covariation is 

[dM k (t),dMj(t)] = (a 2 kju k+:j -i(t) + Ufc+j+i -/3 C fcuk_iUj + i(t) 

(26) - p c ju 3 ^(t)u k+1 {t) + (/3 c ) 2 feiu fe _i(*)u,_i(i)wi(i))dt 

From these covariations we form the covariation matrix S^j(i). Note that this matrix depends upon the 
path of X; the solution to the SDE system (|23|) is conditionally Gaussian given a path of X. 

For the derivation of (|25j) to be rigorous, we must show that X k belongs to the weighted Sobolcv space 
Wq(w,p) for every k <S N. The choice of the weight functions w and p is not unique; see Section [7] We 
will make a choice that is convenient for our purposes but may not be minimal. Let J = 2, p = 1, and 
w = cxp(— cX) for c > 0. Then, p e ~ 1 D e p and w~ 1 p l D i w are bounded for every £ < J and Condition 17.11 in 
Section [7] is satisfied. There exists a unique solution to the stochastic evolution equation for S in this chosen 
space. W J (w,p) is the closure of C^° in the norm || • ||wJ(u>, p v Therefore, to prove A fc belongs to W a J {w,p), 
we show that a sequence /^ € C^j approaches A fc as m — > oo under the norm || • ||vpJ(u> p v 

More precisely, let /*j = X k B m (X) where B m (X) = 1 on [0, if m ], < B Tn < 1 on (K m ,K m + e), and 
-Bm(A) = on [if m + e, oo). Here, (K m ) is a sequence of numbers tending to oo and e > 0. Furthermore, 
B m (X) is smooth. Such a function is 

„ m , I- 00 9m(x)dx 
a m \A) — i -55 , 

Lo°9m{x)dx 

g m {x) = h(x - K m )h(K m +e-x), 

where h(x) = e~* when x > and for x < 0. Then, /^ is a smooth function with compact support. 
We also note that the integral J_ g m {x)dx only depends upon e; it is not affected by K m . To prove that 
\\^ k ~fmW\\w J (w,p) ^ as K m — > oo, one has to show that the following integral tends to zero as K m — > oo 
for each I < J .: 

poo rKm [■K m J r£ 

\ e- cX \D i {X k -fZ l (X)]\ 2 dX= e- cX \D i [X k -X k B m (X)]\ 2 dX+ e - cX \D e [X k - X k B m (X)]\ 2 dX 

JO Jo JK m 

/>oo 

+ / e~ cX \D l [X k -X k B m {X)]\ 2 dX. 

JK m +t 

The first term trivially is zero due to the definition of B m while the third term can be shown to tend to zero 
using integration by parts. Turning to the second integral, we have that 

/ £ e- cX \D e [X k -X k B m {X)]\ 2 dX< [ V cA |D £ A fc | 2 + f ' e- cX \D e [X k B m (X)}\ 2 dX. 
Jk„, Jk„, Jk^ 



The first integral again obviously tends to zero by integration by parts. In order to show that the second 
term tends to zero, one has to demonstrate that the \D e B m (X)\ 2 is uniformly bounded in m on [K m ,K m + e]. 
For I = 0, the result is trivial. For i = 1, we have that 



D 1 B m {\) = Cexp < Cexp 

y{t-y) 



where y = X — K m and y € [0, e]. For £ = 2, we have that 

D 2 S m (A) = Cexpf- 



1 1 

y( e ~y)J V( e ~2/) 2 y 2 

Exponential decay dominates the polynomial growth, so D B m (X) and its derivative approach zero as y 
approaches or e . Since [0, e] is a bounded domain, D 2 B m (X) is therefore also bounded on [0, e]. Therefore, 
since \D l B m {X)\ 2 is uniformly bounded on [K m ,K m + e], the second integral also tends to zero. 

Unless the volatility parameter a = 0, the SDE system ([23]) is not closed^ For computational purposes, 
it must be truncated at some level K. Note that this means one also needs the LLN moments {uk} 2 —^ ■ 
The system (|M)l must also be truncated and can be solved using an Euler scheme. Important details for the 
numerical solution of (f24)) as well as an alternative formulation as a random ODE system are described in 
Gieseckc et al. (2012). 

5.2. Case with no Systematic Risk. When (3 s = 0, there exists a semi-analytic solution. Observe that 
the moment system ([M)l is now a system of ODEs and that (|2"5)l is a linear system with a Gaussian forcing 
term. Therefore, the system of moments v(i) is itself Gaussian. One can directly compute its distribution. 
We rewrite (|25|) concisely as 

(27) dv(t) = A(t)v(t)dt + dM(t), v(0) = v , 

where A : [0,T] ^ rK+i.K+i and v ^ . [o,T] x O ^ R K +\ Now, let * : [0,T] ^ R-K+i.if+i be the 
fundamental solution matrix satisfying 

(28) dtf (*) = A(t)V(t)dt, *(0) = I, 

where / is the identity matrix. If /3 = 0, A is a constant matrix and there is an exponential matrix solution. 
For fj c > 0, one can either solve (f28| numerically using a method such as Rungc Kutta or the Magnus series 
approximation. We assume that Vo is deterministic (if it is a Gaussian random variable, the result is very 
similar). Then, 

(29) v(t) = *(t)v + *(t) f ^- 1 (s)dM(s). 

Jo 

It easily follows that v(t) ~ 7V(*(i)v , E(t)) where 

(30) S(*) =*(*)[ / t *- 1 (.s)S J M(s)(*- 1 (s)) T ^]*(i) T - 

Jo 

Therefore, we have avoided simulation of (|27|) by instead developing a semi-analytic approach for the com- 
putation of the distribution. Also, note that in the case (3 C — f3 s — 0, we have a closed-form formula for the 
approximating loss since ^(t) = e At . 

Figure Q] shows a comparison between the approximate loss distribution according to (|21[) and the actual 
loss distribution for a pool of N = 1, 000 names, for each of several values of the contagion sensitivity j3 . 
Here and in the numerical experiments described below, the actual distribution is estimated by simulating 
the default times in ([T]) using the discretization method detailed in Giesecke et al. (2012). The second-order 
Gaussian approximation (|21[) is a significant improvement over the first-order approximation ([9]) implied by 
the LLN. The latter produces a delta function for this case while the Gaussian approximation is able to 
accurately capture a large portion of the remaining noise in the finite system. 

The linear stochastic evolution equation (|16[) involves a linearization of the jump term /3 dL^ in ([1}. 
It is reasonable to expect that the accuracy of the approximation will decrease as ft increases. This is 
confirmed by the numerical results. For larger pools, the linearization will become more accurate. In the 
case of (5 C = 3, we are taking too small of a pool. 



If a = 0, the LLN reduces to two SDEs and the fluctuation limit also reduces to two SDEs. 
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Figure 1 . Comparison of approximate loss distribution and actual loss distribution in the 
finite system at T = 1 for TV = 1, 000. The parameter case is a = .9, a = 4, A = .2, A = .2, 
and j3 s = 0. 



Another common concern with linearizations is a loss of accuracy over long time horizons. However, we 
do not observe any significant loss of accuracy for the fluctuation limit in our model. 



5.3. Case with Systematic Risk. In the case with systematic risk (i.e., {5 s > 0), the solution to (|25|) is 
conditionally Gaussian given a path of X. There are two approaches to treat the general case. We present 
the most obvious approach first and then a second scheme which achieves lower variance by replacing a 
portion of the simulation with a semi-analytic computation. 



5.3.1. Scheme 1: Direct Simulation. Here we directly discretize and solve the systems (|24[) and (1251) (for 
example, using an Euler method). Alternatively, one could instead numerically solve the random ODE 
formulation of ([24]) given in Gicseckc ct al. (2012) and then solve (|25|) . 

• Simulate paths A 1 , . . . , X M of the systematic risk process on [0, T]. 

• Conditional upon a path X m , the martingale term A4 is Gaussian. 

— First, solve for the LLN moments u m (i) in (|24|) for a path X m and then calculate the conditional 
covariation matrix S^(i). 

— Perform a spectral decomposition of the covariation matrix S^, (t) at a discrete set of times on 
a grid T. 

— For j = 1 , . . . , J, discretize system ([231) on the grid T. At each time ti e T, draw a sample 
dM. m >i(ti) using the spectral decomposition of £^(£) and then solve (|23|) using the Euler 
method. This produces J samples from the conditionally Gaussian solution given the path X m . 



Finally, for each j, compute a sample of the "approximate conditional loss" L 



7n,j,JV 



for a pool 



of size A^ as the difference between the conditional LLN loss LJ" 
the conditionally Gaussian solution scaled by 1/y/N. 
Approximate the unconditional distribution of L^ by -jjj Y] m — 1 X) 7 =i ^i 



1 — UQ l (t) and a sample of 



Let a\ = Var[/ °° f(y)p t (y\X)dy] and a\ = E[Var[/(if r )|X]] where p t (y\X) is the conditional density of the 
approximate loss given X at time t. The total simulation time is Mt\ + MJt 2i where T\ is the time needed 
to simulate the systematic risk process and solve the LLN moment system while t 2 is the time required 
to solve the fluctuation moment system. If we want to estimate the expectation of a function / of the 
approximate loss at some time t, it is straightforward to show using the first order condition for a minimum 
that the optimal allocation of simulation resources M* and J* to minimize the estimator's variance for a 
fixed computational time r is 



M ~ —. — 7; = — ^r T = M n T 



2{o\t x t 2 -o\tI) 



H) 



(31) J* 



1 - M *n 



r 2 A/ * 
The optimal number for J does not depend upon the total computational resources. 

5.3.2. Scheme 2: Conditionally Semi-Analytic Approach. We propose an alternate scheme that replaces some 
of the Monte Carlo simulation of the previous scheme with a semi-analytic calculation. First, we rewrite 
(f2"5")) concisely as 

(32) dv(t) = A(t)w(t)dt + Bv{t)dX t + dM{t), v(0) = v , 

where A(t) is a random, time-dependent matrix and B is a constant matrix. In the calculations that follow, 
we assume that Vo is deterministic; the computations when Vo is a Gaussian random variable are similar. 
The fundamental solution * : [0,T] x i-> M. K+1 ' K+1 satisfies 

(33) d*(t) = A{t)m{t)dt + B^(t)dX t , *(0) = I. 
Then, since [dX t ,dMk(t)} = for any k > 0, 



(34) v(t) = *(t)v + *(t) / ^~\s)dM{s). 

Jo 

(It is easy to show that this is a solution to (|32|) by substituting it back into that SDE.) Unfortunately, even 

when (3 C = 0, bo(x) = 6, and cq(x) = Co so that A is a constant matrix, there is no exponential solution to 

(|33| since A and B do not commute. 

To calculate the approximation (|21j) . one must simulate from the joint law of the paths of X and M. over 

[0, T]. We can do this using the following scheme: 

• Simulate paths X 1 , . . . , X M of the systematic risk process on [0, T]. 

• Conditional upon a path X m , we have a Gaussian solution for the fluctuation moments v m (i). 
Furthermore, we can semi- analytically compute the conditional distribution of v m (i). 

— Given a path X m , solve for the LLN moments u m (t) (which also yield the conditional LLN loss 
L m ) and the fundamental solution \& m (£). 

— Then, form the conditional covariation matrix S^j(i). Finally, compute the conditional covari- 
ance matrix using the closed-form formula ([30]) . This yields Var[z;™(£)]. 

• Approximate the unconditional distribution of Lf by jj X) m =i ^i™> wnere P™ ^ s a Gaussian measure 
with mean L m and variance Var[w™(i)]/iV. 

A skeleton of the systematic risk process X for both schemes discussed above can be generated exactly 
(without discretization bias) using the methods of Beskos & Roberts (2005), Chen & Huang (2012), or 
Giesecke & Smelov (2012). It is also worthwhile to highlight that cither scheme yields an approximation to 
the distribution of the loss L^ for all time horizons t < T (i.e., the "loss surface") and all portfolio sizes 
N G N simultaneously. 

Scheme 2 will have lower variance than Scheme 1. Let f 1 and f 2 be estimators of E[/(L^)] using Schemes 
1 and 2, respectively We have 

M J 



Varf/ 1 ] = Var 



— W f(T m 
JM 






^ + f!)>Var[f]^. 



Var[/ 2 ] = jjo-f since Scheme 2 generates samples from the random variable L f(y)pt(y\X)dy. 
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Figure 2. Comparison of approximate loss distribution for different truncation levels at 
T = .5. Parameter case is a = .9, a = 4, Ao = .2, A = .2, (i c = 1, and (3 s = 1. The systematic 
risk X is an OU process with mean 1, reversion speed 2, and volatility 1. Approximate loss 
distribution computed using the conditionally semi- analytic Scheme 2. The fluctuation 
moments are truncated at level K and the LLN moments are truncated at level 3K. 



We note that there is numerical instability for large /3 for both schemes, especially over long time horizons. 
One must use a small time-step to avoid this instability. The instability is caused by the exponential growth 
terms i(/3 s <7o(At)) 2 fc(fc — l)vk(t) and \{fi s ao(X t )) 2 k(k — l)ufe(i). When one is only interested in calculating 
the LLN approximation ^ via the moment system (|24p. the following transformed moments significantly 
reduce instability by removing the exponential growth term: 



^((3 s ) 2 k(k-l)J <ro(X 8 ) 2 ds 



w k (t) = cxp 



However, when interested in solving (|24[) and (|25[) in conjunction, there is no simple transformation. The 
best approach is to solve ([23]) with a sufficiently small time step such that it is stable and then solve for the 
transformed fluctuation moments Wk(t) = cxp(— i(/3 s ) 2 fc(fc — 1) J Q ao(X s ) 2 ds)vk(t) with a larger time step. 

Figure [2] compares the approximate loss distribution according to (|21| for different truncation levels K 
of the fluctuation moment system (|2"5)l . We use a time step of .005 and produce samples from X using an 
Euler scheme. The approximate loss distribution converges very rapidly in terms of the truncation level. 
This conforms with previous numerical studies of the LLN moment system ([24]) which also demonstrated its 
quick convergence rate; see Giesecke et al. (2012). 

We now study the validity of the second-order approximation pip by comparing it against the actual loss 
distribution from simulating the true finite system. In the following numerical studies, we choose a time step 
of .005 and a truncation level of K = 6 for both schemes described above. A time step of .005 is used for 
simulating the finite system. Samples from X are produced using an Euler scheme. Figure [3] compares the 
approximate loss distribution according to ([21]) with the actual loss distribution for f3 c = and f3 s = 1, for 



35 



30 



25 



20 



15 



10- 

















■ First-Order Approximation 
- Second-Order Approximation - 


: / N =V 

71000 y 


- 


/ N=250 V 


- 


f / \ \ 
,V N=150 v\ 


- 


///■ ^ 




/ft I : 




/ft J : 
// jt 1 : 

// / I .' 

yy / I '■ 




__^-^< y , ■' 


'^V^ir-S^s-^^^^ 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

Loss 



Figure 3. Comparison of approximate loss distribution and actual loss distribution in the 
finite system at T = .5. Parameter case is a — .9, a — 4, Ao = .2, A = .2, j3 c = 0, and /3 s = 1. 
The systematic risk X is an OU process with mean 1, reversion speed 2, and volatility 1. 
Approximate loss distribution computed using the conditionally semi-analytic Scheme 2. 



each of several portfolio sizes. The approximate loss distribution is extremely accurate, even for a very small 
pool with only TV = 250. We also observe that the approximation accurately captures the tails of the actual 
loss distribution. Figure U compares the approximate loss distribution with the actual loss distribution for 
(3 = 1 and j3 s = 1. The approximate loss distribution is again accurate, although not as accurate as when 
f3 c = 0. We also show in both Figures [3] and 0] the first-order LLN approximation ([§]). It is clear that the 
second-order approximation has increased accuracy, especially for smaller portfolios and in the tail of the 
distribution. Finally, Figure [5] shows a comparison for the 95 and 99 percent value at risk (VaR) between the 
actual loss, LLN approximation ([9]), and approximation ([21]) for a pool of N = 1,000. The approximation 
for the VaR based on ([2T]) is significantly more accurate than the LLN approximation for the VaR. 

Appendix [A] provides additional numerical results on the performance of the numerical schemes. In 
particular, we demonstrate the advantages of the second-order approximation over direct simulation of the 
finite system <[T|). For a fixed computational time, the standard error for the second-order approximation is 
several orders of magnitude smaller than the standard error for the finite system simulation. 

5.4. Approximate Loss Process. One of the significant advantages of the fluctuation analysis presented 
here is the dynamic approximation it provides for the entire loss process. An approximation of the loss 
process may be used to estimate the probability that the loss per month never exceeds a certain amount, or 
the density of the hitting time for the loss exceeding a certain level. 

One can generate a sample skeleton of the approximate loss process according to the right hand side of 
([21]) using Scheme 1. Although Scheme 2 gives the loss distribution at all time horizons, it does not directly 
yield a skeleton of the approximate loss process. However, due to the process being conditionally Gaussian, 
one can use a slightly modified form of Scheme 2 to simulate such a skeleton on as fine a time grid as 
desired. For small /3 s , one can expect that Scheme 2 will be more computationally efficient than Scheme 1 
(see Appendix [Al for numerical results). 
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Figure 4. Comparison of approximate loss distribution and actual loss distribution in the 
finite system at T = .5. Parameter case is a — .9, a — 4, Ao = .2, A = .2, (3 C = 1, and j3 s = 1. 
The systematic risk A is an OU process with mean 1, reversion speed 2, and volatility 1. 
Approximate loss distribution computed using the conditionally semi-analytic Scheme 2. 



To simulate skeletons for the moments (f2"5")) of the fluctuation limit, one can use the following method. 
Let v ~ (the general case is very similar). As usual, begin by simulating a path from the systematic risk 
process A. Given a path of A, we then simulate the loss at t = T from the conditionally Gaussian distribution 
given by the formula (|30[) . To simulate at < s < T, we only need the covariancc matrix Cov[vt, v s |Vt]- 
One can continue to simulate at grid points of one's choice. For instance, one can simulate at time si where 
s < si < T given the covariancc matrices Cov[vt, v s |Vt], Cov[vt, v s JVt], and Cov[v s , v Sl \Vt]- These 
covariance matrices are easily computed using the fundamental solution (|33|) : 

(35) E(T 1 ,T 2 ) = Cov[v ri) v„|Vr] = *(Ti)[ H ^ *- 1 (s)S A <( S )(vI/- 1 ( s )) T d s ] *(r 2 ) T . 

Jo 

To simulate v s given v Tl and v T2 where T\ < s < t 2 , we generate a sample from the multivariate Gaussian 
distribution Af(ft, E) where ft = EiEj ^ an d ^ = ^( s > s ) ~ EiE^ ^i- The matrices Si, E 2 , and Z arc 



Si 
E 2 



[S(s,ri) E(s,t 2 )] 

S(ti,ti) E(ti,t 2 ) 
E(ti,t 2 ) E(t 2 ,t 2 ) 



It is worthwhile to note that this scheme requires only a single computation of the fundamental solution 
\E' m (£) for each path X m . Given *I/ m (£) for < t < T, one can generate as many skeletons of the fluctuation 
moments as desired for that particular X m via the aforementioned method. In contrast, Scheme 1 requires 
one to resolve the system of SDEs (j25|) for each new realization. In particular, this suggests that this scheme 
will be more efficient than Scheme 1 for small /3 s . 
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Figure 5. Comparison of approximate and actual value at risks for N = 1, 000. Parameter 
case is a = .9, a = 4, Ao = .2, A = ■2 1 (3 C = 1, and (3 s = 1. The systematic risk X is an 
OU process with mean 1, reversion speed 2, and volatility 1. Approximate loss distribution 
computed using the conditionally semi-analytic Scheme 2. 



Once skeletons from the fluctuation limit have been simulated, skeletons of the approximate loss process 
can be produced by combining these fluctuation paths with the law of large numbers limit process given by 
(f2"4"| and using the approximation formula (|21[) . 

Figure [6] shows the difference between the actual multivariate density and the approximate multivariate 
density for the loss at two horizons. The figure suggests that the second-order approximation is accurate for 
the multivariate distribution of the loss at a collection of times. 



5.5. Nonhomogeneous Portfolio. The moment method can be extended to a nonhomogencous portfolio. 
Define the moments Vk(t, p) = / R+ X k E t (dX, p) and Uk{t, p) = L+ X k fit(dX, p). Then 

"1 



dv k {t,p) = p c ku k _ 1 {s,p)[ / Vl {t,p)dp)dt + 



-cr z k(k- l) + aXk + k/3 L ' I u\{s,p)dp 



v k -i{t,p)dt 



-ka + kf3 s b (X t ) + h(k - 1) (/3 s a (X t )Y 



+ v k (t,p 

+ k/3 s a {X t )v k (t, p)dV t + M k {t, p 
u fc (0,p)= / X k E (dX,p). 



dt - v k+ i(t,p)dt 



The covariation [dM k (t, pi) 7 dMj(t, P2)] for pi, p 2 G V is given by 

[dM k (t, pi), dM 3 {t, p 2 )] = l Pl=P2 (o 2 kju k+j -i(t, pi) + u k+j+ i(t, pi) - /3 c fcii fe _i(t, pi)uj+i(t, pi) 
- f3 c ju j ^ 1 (t, pi)w fe+ i(t, pi) + (I3 c ) 2 kiu k -i{t, pi)wj_i(t, pi)«i(*, p\))dt. 





Figure 6. On the left, the figure shows the difference between the actual multivariate 
density and the approximate multivariate density for the loss at T = .25 and T = .5 in a 
pool of N = 1, 000 names. On the right, the approximate density is shown. Parameter case 
is a = .9, a = 4, A = .2, A = .2,/3 c = 1, and (3 s = .4. The systematic risk X is an OU 
process with mean 1, reversion speed 2, and volatility 1. 



The LLN moments are 

du k (t,p) = u k (t,p) 

+ Uk-i{t,p) 



-ak + (3 s b Q (X t )k + -(/3 s ) 2 a 2 Q (X t )k(k - 1) 



0.5cr 2 fc(fc- l)+aXk + /3 c k / ui(*,p)dp 

v 



dt- u k+1 (t,p)dt 

dt + p s <7 (X t )ku k {t,p)dV t 



u fc (0, p) = / A> (dA, p). 
Jo 



For numerical implementation, one must discretize the parameter space V . As we allow more parameters 
to vary, V can become high dimensional. This can become computationally expensive. 

6. Proof of Theorem 14.11 

In this section we prove the fluctuation Theorem 14. II The methodology of the proof goes as follows. After 
some preliminary computations, we obtain a convenient formulation of the equation that {5^} satisfies, see 
(|4ip . Some terms in this equation will vanish in the limit as N — > oo; this is Lemma |6. II Based on tightness 
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of the involved processes (the topic of Section [5]) and continuity properties of the operators involved, we can 
then pass to the limit and thus identify the limiting equation. The limit process satisfies in the weak form 
another stochastic evolution, which has a unique solution (proven in Section [5]). The difficulty is to identify 
a rich enough space where tightness and uniqueness can be simultaneously proven. This is not trivial in 
our case, because the coefficients are not bounded and the equation degenerates (coefficients of the highest 
derivatives are not bounded away from zero) . It turns out that the appopriate space in which one can prove 
both tightness and uniqueness is a weighted Sobolev space, introduced in Purtukhia (1984) and further 
generalized in Gyongy & Krylov (1990) to study stochastic partial differential equations with unbounded 
coefficients, and are reviewed in Section [7J 



An application of Ito's formula shows that for / <E C£(V), 



(36) 
where 

(37) 



+ f A N [f]{X s ,^)ds+ f(c*'f,tf) E dV a + (f,M?) 



ds 



N 



J, 



f 



n=l 
/ \ 1 V^ I t I N.n' \N,n' | PN,n\ f ( N.n' \N.n'\ { ,,N,n' 1 t(„N,n \N,n\ 



for all s> 0, N £N and n£ {1,2,..., N}, 

N 

(38) 



and 
(39) 



n—l " n— 1 ^ 

M?> n = (1 - Mf >») - 



\^ n M?' n dr. 



By subtracting ((5J from (|36|) we find that S^ satisfies the equation 

</, ~f ) = (/, ~ N ) + jf* (A/ + £f /, ~f ) ds + £ [<Q, H.) (C 2 f, Ef ) + (C, Hf ) (£ 2 /, M f >] ^ 

/NA N [f](X s ,^)ds + J (c?'f,E?)dV a + (f,VNMl 



In order to simplify the calculations later on, for each / £ C$(V), t > 0, N £ N and n £ {1,2,..., AT}, we 
define the quantity 



- r f ^ dof PjV,n V^ ®f 



Jin(t) 



EVJ ft .N,m\ 



AT ^ 0A 

m— 1 



M 



N.m 



f(p N ' n ,\?' n ) =ffi, n (c 2 f,n?) E -m?). 



A simple computation shows that 
(40) 



•^N,nW \T<^N,nW 



< 






a 2 / 



9A 2 



c 



We can rewrite the stochastic evolution equation that E N satisfies as follows: 
(/, Sf ) = (/, ~ W ) + J* (A/ + £f /, ~f ) ds + jf [(Q, jS.) (C 2 f, Sf ) + (Q, ~f ) <£ 2 /, M f >] d s 

(41) + J*(c?'f,E?)dV. + (f,VNM?)+IlZ 
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where 



/, Vnm? 



-(is/: 



&N.: 



A * — (P» > M s dW s +^1^ 



H t,o - 



N 



JnJ s )^ 



N,n 



</NA N [f](X s ,^)ds + f^ [ VNB N ' n [f](s)dM^ n 



B 



N,: 



1 s, 



[/](») = ^»-Tf^»- 



AT 



The term i?^ turns out to vanish in the limit as the following lemma shows. 
Lemma 6.1. For any t £ [0,T] and any f £ C%(V), there is a constant Cq, independent ofn,N such that 



and 



Moreover, we have that 



3 



sup E 

0<t<T 



N- 



A"[f\(X r ,tf) 



dr 



< 



Co 



V / ^NB N > n [f](s)dAf s N > n 



lim E sup K o | 2 =0. 

JV-s-oo <t<T 



2V 



- AT2 



Lastly, regarding the conditional (on V t ) covariation of the martingale term ( /, yNAi^) we have 

ftiAt 2 



Cov 



f,VNMZ),<g,VNMi 



<N 



V tl N 



V, 



(C B (f,g),ij?) + (&(/, g),tf) + 

(42) - <£ 7 <7, /xf ) <£ 2 /, /if) - <£ 7 /, Mf ) (^25, ^ >] ds 

We continue with the proof of the theorem and defer the proof of Lemma [B~T1 to the end of this section. Rel- 
ative compactness of the sequences < S t '",£ e [0, T] > and < y/~NM t ' n ,t £ [0,T]> iriD -t_j+ 2 ), JO,' 

follows by Lemmas l8.7l and f8.8l Relative compactness of the sequence {fj,^ , t £ [0, T]~j N „ in -D_e[0, T] follows 
by Lemma 7.1 in Giesecke et al. (2012). These imply that the sequence 



tlVt 2 



{(<•"' 



M N ,E N ,V\ ,N£N\ 



is relatively compact in D ExW - i j+2) {wp)xW -(j+2 ){wp)xR [Q,T}. Let us denote by 

{(ji,M,E,V)} 

a limit point of this sequence. We mention here that this is a limit in distribution, so the limit point may not 
be defined on the same probability space as the prelimit sequence, but nevertheless V (and thus X), have 
the same distribution both in the limit and in the prelimit. Then, by (|41|). Lemma 16.11 and the continuity 
of the operators Q and £4 we get that { (fi, M., S, V) } will satisfy, due to Theorem 5.5 in Kurtz & Protter 
(1996), the stochastic evolution equation (fTB"]) . By Lemma IB~T1 the conditional covariation of the distribution 
valued martingale Ai is given by (|18|) . Uniqueness follows by Theorem 19.71 with J = 2. This concludes the 
proof of the theorem. 

We conclude this section with the proof of Lemma 16.11 



Proof of Lemma \6.1\ We start by recalling that Lemma 3.4 in Giesecke et al. (2013) states that for each 
p > 1 and T > 0, 

N 



(43) 



<t<T JV ' L 



0<t<T 

jYt:i 



<C 
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for some constant C > that is independent of TV, n. Along with (|4"0")) , this implies that 



E 



A»[f\(X r ,tf) 



dr 



E 



< K z 



^o „ =1 iV n=1 

1 W 
wE A ^ 



f/v 



a 2 / 



< 



c, 



iV 



/ 



for some nonncgative constant Co that is independent of N, n. Then, the first statement of the lemma 
follows. The second statement follows similarly. In particular, using the martingale property, pt)|) and (|4"3")l 
we immediately obtain 



E 



' N 

E 

,n=l 



-, 2 



7VB iV ^[/](s)dAA s f 



N.n 



Co 

~ N 2 



for some constant Cq > 0. The statement regarding R^ follows directly by the previous computations and 



Doob's inequality The statement involving the conditional covariation of the martingale term ( /, 
follows directly by the definition of Ai N and the fact that jumps do not occur simultaneously. 

7. The appropriate weighted Sobolev space 



NM¥ 



D 



In this section we discuss the Sobolev space in which the fluctuations theorem is stated. Similar weighted 
Sobolev spaces were introduced in Purtukhia (1984) and further generalized in Gyongy & Krylov (1990) to 
study stochastic partial differential equations with unbounded coefficients. These weighted spaces turn out 
to be convenient and well adapted to our case. See Fernandez & Meleard (1997) and Kurtz & Xiong (2004) 
for the use of weighted Sobolev spaces in fluctuation analyses of other interacting particle systems. 

Let w and p be smooth functions in V with w > 0. Let J be an integer and define by Wq(w : p) the 
weighted Sobolev space which is the closure of C X> ('P) in the norm 

1/2 



l/ll 



Wg(w,p) 



,fe<J 



w'(p) p«(p)D«f(p) dp 



< oo 



pS"P 



For the weight functions w and p we impose the following condition, which guarantees that Wq (w, p) will 
be a Hilbcrt space for every J > 2, see Gyongy & Krylov (1990). 

Condition 7.1. For every k < J, the functions p k ~ 1 D k p and w^ 1 p k D k w are bounded. 



The inner product in this space is 

\fi9)wJ{w, P ) 



E 

k<J 



w 2 (p)p 2k ip)D k f(p)D k g(p)dp 



i>ev 



Moreover, we denote by W (w,p), the dual space to Wq(w,p) equipped with the norm 



l/ll 



"n" 



sup 



<f,9> 



l/ll 



(w.p) "T ||_|| 

g£W 6 '(w,p) llffllvi^Kp) 

For notational convenience we will sometimes write ||/||j and ||/||_/ in the place of II /ll vk- 7 (u; p) 



and 



W ~ J (,w,p) 



respectively, if no confusion arises. 



The coefficients of our operators have linear growth in the first order terms and quadratic growth in 
the second order terms. For example when J = 2, we easily see that the choices p(p) = \J\ + |p| 2 and 

w(p) = (l + |p| 2 ) with (3 < —3 satisfy the assumptions of Condition 17.11 

In the next lemma we see the necessity for introducing weighted Sobolev spaces. Notice that if the weights 
w, p are chosen as above, then the condition on p of Lemma 17.21 is equivalent to assuming that /i££ has 
finite moments up to order 2/3. 
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Lemma 7.2. Let Condition \7.1\ hold and fix x £ K and p <G E such that Jp (u> 2 (p)p 2 (p)) p{dp) < oo. 
Then the operator Q x ^ is a linear map from W Q (w, p) into Wq (w, p) and for all f £ W (w, p), there 
exists a constant C > such that 



2 
lj+2 



||^/|| 2 7 < C (l+ \b (x)\ 2 + \a (x)\ 4 + (&, p) 2 ) ||/|| 
Moreover, the operator C\ is a linear map from W (w, p) into Wq (w, p) and for all f £ W (w, p) 

WClffj < C\a Q (x)\ 2 \\f\\ 2 J+1 < C\a a (x)\ 2 \\ff J+2 
Proof. Recall the definition of Q x ,p,f and examine term by term. For the first term, Condition l7. II gives us 



i^i/ii 2 / = E/voV fc (p)(^i/(p)) 2 rfp 

k=o Jv 
< C J2 l,w 2 {p)p 2k {p)(D k f{p)) 2 dp 



fc=0 

< C\\f\ 2 



■7+2 

P 
,7+2 



For the second term again Condition 17.11 gives us 



I^/II* = Y,(w' 2 {p)p 2k {p){D k Llf{p)) 2 dp 



J+2 . 

< C(|& (s)| 2 + MaO| 4 )£ / w 2 (p)p 2k (p){D k f(p)) 2 dp 

k=o Jli 



< C(\b (x)\ 2 + \a (xt)\\f\\ 2 j, 
For the third term we have 



J . 
Qi,n)£2f\\ 2 j = (Qi,M> 2 E / w 2 (p)p 2k (p){D k C 2 f( P )) 2 dp 

k=v Jv 



< c(Q llM r ii/iis+2 

For the fourth term we have 



C2f,li)Qifj < C(C 2 f,p) 2 J2 f w 2 (p)p 2k (p)(l + X 2 )dp 

k =o Jp 

< c(JjDf(p)\n(dp)) Y^J^w 2 {p) P 2k {p){l + \p\ 2 )dp 

< f w 2 (p) P 2 (p)\Df(p)\ 2 dp f ( w 2 ( P y( P ))-V(#)E i w 2 {p) P 2k (p){i + \p\ 2 )dp 

JP JP ,.,_n JP 



k=0 

< C\\f\\ 2 j 



(44) < C||/|| 2 J+2 

this last computation shows the need for the weighted norm. The statement for C% follows analogously. This 
completes the proof of the lemma. □ 

8. Tightness and continuity properties of the limiting process 

We first recall three key preliminary results from Giesccke et al. (2013) and Giesccke et al. (2012) that 
will be useful in the sequel. 

1!) 



Lemma 8.1. [Lemma 3.4 in Giesecke et al. (2013).] For each p > 1 and T > ; there is a constant C > 0, 
independent of N, n such that 

N 



sup ^£E[|AfT 



0<i<T 

JViEN 



71=1 



<c. 



Lemma 8.2. [Lemma 8.2 in Giesecke et al. (2012).] Let W* be a reference Brownian motion. For each 
p = (p, A ) € V where p = (a, X, a, (3° , (3 s ), there is a unique pair {(Q(t), Xt(p)) '■ t <E [0, T]} taking values in 
R + x R + such that 

(45) Q(t) =y*j3 c E (VjW JA t *(p)cxp -y A.:(p)d S 
and 

(46) A*(p) = Ao - a / (A:(p) - A)ds + cr / v^»™C + / Q(s)ds + /3 s / X* s (p)dX s 



v{dp). 



Lemma 8.3. [Lemma 8.4 in Giesecke et al. (2012).] For all A € <3S(v) and B € 38(M. + ), \i satisfying f3|) 
is given by 



fjL t {A xB)= X A(p)E (Vt ,p) 



Xs(A t *(p))exp 



X* s (p)ds 



v(dp), 



where Ajf(p) is defined via J^5[j-|7^). Therefore, for any f € C(V), 



(f,lH) E = / E (v t .p) 
J-p 



/(pA*(p))exp 



X* s (p)ds 



Let us now set 



Notice that 



G(t) d = f jj c (l E (Vl p) 



exp 



A:( P )d S 



i/(dp). 



f(dp). 



G(<) = Q(i). 



We prove tightness based on a coupling argument. The role of the coupled intensities will be played by 

c tN 

N.n-^t 



X t '"' , defined as follows. Define A t '"' to be the solution to ([TJ when the loss process term /3^ n L^ in the 
equation for the intensities has been replaced by G{t). Let us also denote by T N,n,G the corresponding default 



-N,G 



time and by L t ' the corresponding loss process. It is easy to see that conditional on the systematic process 
X, the \ N - n - G for n = 1,...,N are independent. The related empirical distribution will be accordingly 
denoted by \i t ' . 

Let k > and define the stopping time 

(47) 



C 1 N 

I, 71— 1 



+ |Af ' n ' G \ p 



> K P 



Notice that the estimate in Lemma 18.11 implies that for T > and p > 1 



(48) 



lim sup 



• N£N 



N. 



,„ < r) = o 



fr, JV,1 



Thus, it is enough to prove tightness for {S tA '£ , i S [0,T]}jvgn- 
The following lemma provides a key estimate for the tightness proof. 

Lemma 8.4. Let J = 2. For every T > 0, i/iere is a constant C independent of N such that 

2 



sup sup E ( N 

N£N0<t<T 



^tA0jv,„, 



N,G 



W^\w,p) 



<c. 
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Proof. For notational convenience we shall write 9 n in place of 8n,u,k- Let us denote by rj^ = /j,^ — /j, t ' 
and let / £ C^(P). For notational convenience we define the operator 

(&,*,m/)(p) - (A/)(p) + (£|/)(p) + (Qi, i/) (£ 2 /)(p) + (£2/, m) Ci(p) 

After some term rearrangement, Ito's formula gives us, via the representations of Lemmas 18. 2118.31 



<rf> 



(Gx.,tf-°,n»f>1?)<l» + J (C 2 f,^' G )(Qi,^' G -ps)ds + j (cf'f,r^)dV a 

1 N ft N ft ft 

^E f WW?" ~ JrJt, /V(Af- G )^f- G 
n=l^° n=l"' 

^ E f (°s„yJ>&c*fW) - a N , n ^^C,f(X^- G ?j dW? 



N 

Using the bounds of Lemma 16. 1[ Ito's formula and the fact that common jumps occur with probability 
zero, we have 



E(/.^ n > 



*M ) (Sx.,tf-°, lt ?f>n?) + \(£*'f>rf 



X{6 n >s} 



ds 



+e£ 2 (/, tf) (£ 2 f, ^ G ) (Qun?' G - fi s ) X{ e n > s] d S + O (1) 



where the O (jj?) term originates from the martingales and from Lemma 16.11 Next we apply Young's 
inequality to get 



Jo 



t(f>n?){Sx.,tf<°, l >»f rf) + 



Vt m fX 



2(f,r,?) 2 + 1 -(C 2 f,^> G ) 2 (Q 1 ,^> G -fi 8 ) 2 



X{e n > s }ds 



X{e n > s }ds + O [ — 



E 



Hf>n?)(0 x .rf«rff,n?) + \(c?'f,rfi)\ 2 + 2(f t n?) i 

-L / „ at n\ 2 / ~ n n — \ 2 , ^ / t 



X{e„> s } 



rf.s 



+E / - <£ 2 /, ^^ (fir, ^ G - /*.)' X{ e n >s}ds + o[- 







AT 



Let us next bound the term on the last line of the previous equation. It is easy to see that 



n d = f ^V<Qi,^ G -^) = E(/0 



N,n,G M N,n,G 



-{Qufi.)) 



is a discrete time martingale with respect to P(-|V). Thus, by Theorem 3.2 of Burkholder (1973), we have 
that 



(49) 



E 



r. 



v 



< CE 



\S(Y S 



V 



where, maintaining the notation of Burkholder (1973) 

S(Y S ) = 



N 

E K 



\ n=l 



Ag M« - (fil, £is 



Therefore, recalling the bound from Lemma 18.1 1 we obtain 



^\Y S \ 2 


= 1e 

N 


E 


'\Ys\ 2 


V 






< ^E 

~ N 


E 


)S{Y S )\ 2 


vj 
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^E^[^ n Af- G Mf^-( 2l ,/l s 



(50) 

This bound implies that 
1 



C 



E 



<£ 2 /,^'T<Qi,/^' G -/i.> X { e n > s} 



d s <-C(t,K, K )\\f\\^ M 



Hence, if {f a }a>i is a complete orthonormal basis for Wq (w, p) we obtain that 



E(/a,7^ n ) < / E 

Jo 
Summing over a > 1 we then obtain 

el 



^{U,rff)(Q Xttll? ,G^f a ,rff) + \(cf°f a ,r,?)\ 2 + 2(f a ,r^) i 



1 2 

X{e„> s }rfs+ — C ||/o|lw J(t(,,p) 



MVtke 



< / E 



2 \^' a x a , /i f'°,„f^)_ J + 



r *,X s N 



+2 k 



ATI 



X{e„>s}rfs+ ^C- 



Then as in Lemma 19.61 we get that 

E 



v N - G 1> 



2(n s ,G~ n,g llN V, 



r*' X 'n N 



K\K*»S-j 



X{e n >s} < CE||?7 S \\_jX{9 n >s} 

1 



ds + 



N 



®\\r& J_j<C 
Finally, an application of Gronwall's lemma concludes the proof. 

The following lemma provides a uniform bound for the fluctuation process. 



□ 



Lemma 8.5. Let J = 2. For every T > 0, there is a constant C independent of N such that 

2 
sup sup E 



N£N0<t<T 

Proof. Clearly Sf ■= VN(jj^ - (X?' G ) + VN(jj^' G - p t ) 

N 



<c. 

Wf J (w,p)/ 

* 1 ,/atN,G 



Nrjt + v Nr/ t ' . Notice now that by Lemmas 



and 



N \f' , h / = 2_ [f( Xt ' ) M 



N,n,G 



(f,fit) 



is a P(-|Vt) discrete time martingale, which in turn implies via Theorem 3.2 of Burkholdcr (1973) (similarly 
to the proof of Lemma I8.4|) that 

*{frf- a y<j;c\\n w , {UtP) . 

This estimate, Lemma T8.4I and Parseval's identity yield the result. □ 



Next we discuss relative compactness for {\Z~NjCi^ ,t £ [0,T]}jvsn- In particular, we have the following 
lemma. 

Lemma 8.6. Let T > and J = 2. The process {vNAif , t £ [0, T]}jvgn is a W Q ~' (w, p) — valued martingale 
such that 

2 



(51) 



sup E 

JVeN 



sup 

0<t<T 



NM 



N 



Furthermore, it is relatively compact in D -(.1+2) , . 



W- J (w,p) 

[0,21. 



<C<oo. 



Proof. Let {f a } a >i be a complete orthonormal basis for W (w,p). Clearly, it is enough to prove tightness 



for {VNM^. g p ,t£ [0,T]}7v eN . Let us define 



(52) Ff [/] = (£ 5 /, „?) + {C 6 f, tf) + (C 2 f, p?Y <Q 2 , tf) 2 (£ r f, tf) (C 2 f, tf) 

Then, an application of Lemma 18.11 shows that there exists a constant C such that 
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(53) 

Therefore, we obtain 



r s [/k{6 



t > s > 



<is 



& N 



<C(t,K, K )\\ff wiM (t-r) 



E 



NM 



N 



NM 



N 
r/\8* 



&N 



W- J (w,p) 



\-? 



< 



E> 



f a ,VNM", 



NM 



N 



&N 



\-? 






(54) 



Similarly, we can also show that 



< C{T,K,k) 



zLWfoWw,! 



O.p) 



a>l 



(t-r) 



3 



sup 
o<t<r 



NM 



N 
tA6* 



";; 



>(™,p) 



< C(T,K,K)Y,\\fa\\w. ( 



O.p) 



(55) 



supl 

AT 



< 



'/ 



Let us define Bs = {/ £ Wq J (w, p) : WfWw'tw P ) — $}• Using the previous bound, we obtain that for 
/ £ Bs, and for every T > 0, r) > 0, there exists a constant tt^T, n, 6, K, k) such that 

sup (f,VNM& g , Y >m(T, V ,S,K, K ) 

_0<t<T X *>,*,"/ 

, t £ [0, T] > satisfies the compact 

containment condition in Wq (w,p), which together with (|54|) gives relative compactness of {y/NM^ ,t £ 
[0,T]}at gN in D w -j, p) [0,T] (Theorem 8.6 of Chapter 3 of Ethier & Kurtz (1986)). 

The uniform bound of the lemma follows by the previous computation and by the fact that 9 P N 
as n — > oo, see (|4"8l . 



Since Bs is a compact subset of Wq(w,p) we obtain that <^ y/NM^.gp 



Next, noticing that 
follows. 



] w r : 



oo 



l{wp) .- - n \\ w - J (w. P y the claimed tightness of {VNM? ,t £ [0,T]}Ar eN 

D 



<C| 



Regarding the convergence of the martingale VNM N we have the following lemma. 

Lemma 8.7. LetJ = 2. The sequence <\/NMf ,t £ [0,T}> is relatively compact in D w -( j +2), , [0, T]. 

Moreover, it converges towards a distribution valued martingale {Mt,t £ [0,T]| with conditional (on the 
a— algebra V) covariance function, defined, for f,g £ Wq (w,p), by U8]). The martingale {Mt,t £ [0, T]} 
is conditionally on the a '— algebra V, Gaussian. 

Proof. Relative compactness follows by Lemma [8.61 Due to continuous dependence of (|4"2")l on p, N and on 
the weak convergence of jjl n — ¥ ft. by Giesecke et al. (2012), we obtain that any limit point of \/NM^ 



as N — > oo, M, will satisfy (fT8[) . Conditionally on V, the limiting M is a continuous square intcgrable 
martingale and its predictable variation is deterministic. Thus, by Theorem 7.1.4 in Ethier & Kurtz (1986), 
it is conditionally Gaussian. This concludes the proof. □ 



Next we discuss relative compactness of the process I B 4 '", t £ [0, T] > 



NEN 



r~N,; 



Lemma 8.8. Let T > 0. The process {^ t '" ,t £ [0,T]}jv e N is relatively compact in D _(j +2 ), . [0,T]. 
Proof. Let {f a } a >i be a complete orthonormal basis for Wq (w,p). It is enough to prove tightness for 

By (|4"Tj) we have 
(56) (/„, Hf ) = (f a , Sf ) + /' {Gx^^sfa, Sf ) ds + /' (£f ■ /„, Sf ) dF s + (/„, VNM&) + < r . 
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Next we notice that 

E E 



■z;N\ 2 



(Qx s ,p,„^fa,£*) X{e* n >e}ds\& 



&N 



< 



E E 



I^JVII 2 



ST A 7 



Gx.,p.,^fa\\ w j (w>p) ||3 8 \\w - J (w,p) X ^ N n K > s }ds\^ r 



< C £e f \\f a \\\ > J+2(w!p) \\E?\\ 2 , X{e » N ^> 8} ds\J? r N 



< CE 



I ■— -iV II 2 



&N 



E s \\w ->(w,p)X{8 P N , n , K >s} ds \^r 



where the second inequality is due to Lemma 17.21 Similarly we also obtain that 



E E 

a>\ 



Ci'fa,2?) X K _>s}ds\^ 



&N 



< CE 



I -AM 



\W- J (w, P ) M»N 



*{6 P N n . K >s} ds \^r 



N 



The last estimates, the uniform bound from Lemma [8. 5 1 Lemma [8.61 for y/NJ\4^ and Lemma [6.11 for the 
remainder term Rf r imply the statement of the lemma. We follow the same steps as in the proof of Lemma 
18.61 and thus the details are omitted. □ 

We end this section by proving a continuity result. 

Lemma 8.9. Any limit point of {E^,t £ [0,T]}jveN and {y/NN^ ,t G [0, T]}jveN is continuous, i.e., it 
changes values only in C w -(j+2), J0,T]. 

Proof. In order to prove that any limit point of {S^,£ <E [0, T]} „„ changes values only in C w -(j+2), JO, T], 



it is sufficient to show that 



lim E 

Af->oo 



II— N —AT ,, 

sup p t _ "4-llM/-< J+2 »r 

t<T 



I Wp-" 1 -" (w,p) J " 

Let {f a }a>i be a complete orthonormal basis for Wq +2 (w, p). Then, by definition, we have 

(/ 0! Sf -Sf_) =VN[(f a ,^- Pt ) - (f a , ^.- fit-)] 

= s/N [(f a ,rf - nf_) - (f a ,p t - fit.)] 

ffcjt)- (fa, fit ~ fit-) 



^>N,nV-) jyJN,nV-) 



+ -^Jfrjt) -VN[(f a ,ik - fl t -)} 



Clearly, 1 1— V m is continuous, so we only need to consider the first two terms. The first term is bounded by 
see (|40|) . whereas for the second we immediately get 



a 2 la 

'5x T 



Jfrjt) <K[\\f a \\ + \\f' a \\] 



So, we overall get that 

(57) 

This implies that 

E 



E sup \(f a ,EF-3*L)\<C 
te[o,T] 



1 



1 



7V3/2 Jfi 



\fa\\\Y^ +2 (w,p) 



II— N — N || 
«j, " Uw o (w,p) 



< 



E- 



a>i 

<C 



t<T 



iV 3 / 2 v/F 



/ , ll/ttllli 



W w+ ^(^,p) 
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which, after taking the limit as TV — > oo, concludes the proof of continuity of the limit point trajectories of 

{~?,te[0,T]} Nm . 

Next we consider continuity of the trajectories of the limit points of {y/NN^ ,t G [0,T]}at g n. It follows 
directly by ([56]) and Lemma \6. II that Sf and yfNM^ 1 have the same discontinuities. Thus we obtain 



lim E 

JV->-oo 



sup 

t<T 



NM t N - \fNM?_ 



i In 



>{w,p) 



< 



lim Ve 



N^roo 



sup {f a ,VNM t 

t<T 



N 



^ J™.,4fEii^n^ 



N^OC, y/N 



(w,p) 



NMfL 



which concludes the proof of the lemma. 



9. Uniqueness 



□ 



In this section we prove uniqueness of the stochastic evolution equation (|16|) . Let Ej, S^ be two solutions 
of (fTo]) and let us define $t = Sj — Sj . Then, $ t will satisfy the stochastic evolution equation 

(/,$*}=/ (Cif + /%'f,*,)ds+f [(Q,ji s )(£2f,$s) + (Q,® s )(C2f,Ji s )]d S +[ (c$'f,*.)dV a 

(58) = /" (Qx.,nJ,^.)da+f (c%°f,$ a )dV a a.s. 

Notice that this is a linear equation. In order to prove uniqueness, it is enough to show that 



(59) 



E||$t| 







Let {f a } be an orthonormal basis for W (w, p). By Ito's formula we get that 



K/-,*t>r 



2 (/„, $ fl ) (A/a + £3 ' /a, $s ) + (Q, M.) (£* fa, *«) + (Q, *s> (A/a, Ms 



^/a,^ 



ds- 



(60) 



'ds + J 2(f a ,$ a )(c?'f a ,*,)dV. 
j 2(f a ,$ s ){gx a , fi J a ,$ s )ds+ J |(£f s ./a,$ s )| 2 d s + J 2{f at $.)(C?' /*,$.} dV. 



Then, summing over a and taking expected value (due to Lemma 17. 2[ the expected value of the stochastic 
integral is zero) we have 

ft r-t 2 

(61) 



E||$ f || 2 _ J= E f 2($ s ,g* s $ a )_jds + E j 
Jo Jo 

Hence, if we prove that there is a constant C > such that 



Cl 1 *'*, 



da 



(62) 



2{<j>,Q*<t>)_ J + \\cr4>\\ 2 _ J <c\\<t>\t J 



then we can conclude by Gronwall inequality that <&t = a.s. 
Let us recall now that 



(&,m/)(p) - (A/)(p) + (A/)(p) + (Qi, /i) (A/)(p) + (A/,m> Si(p) 



and set 



(# iM /)(p) = (A/)(p) + (£f/)(p) + (Qi,/x) (A/)(p) 
(5* iM /)(p) = (A/,m)Qi(p) 

and (Glf)(p) = (Qk^JM for i=l,2. _ 

Moreover, for simplicity in presentation and without loss of generality, wc shall consider from now on 
only the case of a homogeneous portfolio, i.e, set v = Sp . This is done without loss of generality, since the 
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operators £,/ involve differentiation only with respect to A. The proof for the general heterogeneous case, 
is identical with heavier notation. 

Before start examining term by term, we gather some straightforward results in the following lemma. 

Lemma 9.1. Let ijj £ C^°(IR+) and J > 1 and assume Condition \ 7. 1\ Then, up to a multiplicative constant 



for the term O I \\tp\\j ) , that may be different from line to line, we have 
J „ 

J2 w 2 (x)p 2k (x)\^ k+1 Hx)^ k \x)dx = o(\\ip\\yj 

j^j w 2 (x)p 2k {x)^ k +'Hx)^ k Hx)dx = o(\\n 2 j) 

jhf w 2 (x) P 2k (x)Xi^(x)^ k Hx)dx = -f w 2 (x) P 2k (x)x i^\x) 2 dx + o(\m 2 j) 

^2 I w 2 {X) P 2k {x)x 2 ^ k+2 \x)^ k \x)dx = -f w 2 {x) P 2k {x) x^ k+1 \x) 2 dx + o(u\\ 2 j 



k=0 JK + 



Proof. It follows directly by integration by parts using Condition 17.11 and the assumption that ip and its 
derivatives arc compactly supported. □ 



Then we bound each term on the right hand side of (|6ip. First, we notice that, by Riesz representation 
theorem, for <f> £ W ~ J there exists a unique ip = F(<p) £ W J such that 

(<f>,f) = (i>,f)j 

By a density argument we may assume that ip = F(<f>) £ Wq and obtain 

(63) (<p,g*<f>) _j = (^,g*<p) = (G4>,cf>) = (gii>,i>)j. 

which is true since by Lemma I7T21 Qip £ Wq . 

Lemma 9.2. For (p £ W ~ J such that ip = F{<j>) £ Wq +2 we have 

{<P^*<P)_j = (l + |6o(x)| + |ao(.T)| 2 + (Q 1 , Ai ))o(||0|| 2 _ J )-i|/3 s a o (x)| 2 ^ / w 2 (X)p 2k (X) |AV^ +1) (A) 
-r'E/ w 2 (X)p 2k (X)X^(X) 2 dX-J2[ w 2 (X) P 2k (X)X ^W(A) 



dX 



(64) 



fe=0" K + 



dX 



Since ip £ W Q , the integrals on the right hand side are well defined and bounded. 
Proof. We recall that 

(65) (<p,g 1 '* ( p)_ J = (g 1 ip,ip) J . 

and that 

(66) (SVXp) = (AVO(p) + (£f V0(p) + (fii, n) (£ 2 V0(p) 

Hence we bound each of the three terms separately. Using the statements of Lemma [97T1 we have for C\ip 



k=o " R + 



W w 2 (X)p 2k (X)i^(X)(C 1 ^)) {k) dX 



\a 2 Y, I w 2 {X)p 2k {X)^ k HX) [x^ 2 \X)) ik) dX-aY: f w 2 (X) p 2k (X)iP^ (X) (x^{X)) W 



fc=0" K + 



dX 



fc=0" M + 



2(i 



+ a\^2 / w 2 (\)p 2k (\)^(\)^ k+ V{\)d\-J2 / w 2 (X)p 2k (X)^(X)(X^{X)) {k) dX 

= l^lt / w 2 (X)p 2k (X)^(X) ( K ^ k+1 \X) + X^ k+2 \X)) rfA 

./ . 

-J2 a w2 ^P 2k M^* 5 ( A ) ( K ^ (fe) ( A ) + A ^> (fe+1) (A)) rfA 

fc=0 ^ R + 

+ E aA / w 2 (A)p 2fe (A)^W(A)^ (fe+1) (A)dA-£] / w 2 (A)p 2fe (A)V (fe) (A) (^^(A) + A^ (fe) (A)) dX 
k=o ^ R + fe=0 ^ R + 

= 0(U\\ 2 j) -\fjl f w 2 (X) P 2k (X)X^HX)\ 2 dX-Y / f w 2 (X)p 2k (X)x\^ k \X) " 



(67) 






fc=0 

J 



Ofll^ll-./) -|- 2 E / ^ 2 (AV 2fc (A)AU fc+1 >(A)| 2 dA-£ f w 2 (X)p 2k (X)Xy k \X) 



dX 



dX 



Similarly, for CZip we have 



fc=0" K + 



ik) 



J2 /_ ^ 2 (A)p 2fc (A)V (fc) (A)(/:^(A)) (fc) dA = 

= /3 s &o(*)]T / «; 2 (A)p 2fc (A)^W(A)(AV^)(A)) W dA 

fe=0 JR+ 

+ ^|/3^o(x)| 2 X: / - 2 (A)p 2fc (A)^(A)(AV 2 )(A)) (fc) dA 

= (3 s b (x)J2 / w 2 (A)p 2fe (A)^( fe )(A)(fc^ (fe) (A) + A^ (fe+1) (A))dA 



fc=0" M + 



A:! 



+ ^o(*)| 2 £ j ( » 2 (A)p 2t (A)^ fc )(A)(^^(A)+2a^n w ^A)^(A)^p^ 



rfA 



dA 



= (b Q (x) + \p s a a (x)\ 2 )o(u\\ 2 r )--\t3 s a (x)\ 2 Y, / ™ 2 (A)p 2fc (A) |A^ fe+1 >(A) 
= (b (x) + \p s a (x)\ 2 )o(uf_ J \ -i\^a (x)\ 2 J2 / ^ 2 (A)p 2fc (A) A^ +1 >(A) 



fe=0 ,/K + 



(68) 



Lastly, for the third term, (Qi,p) (C2tp)(p), we have 



£ 

fc=0-' K + 



(Qi,l*)J2 ^ 2 (A)p 2fe (A)^ (fc) (A)(/:2^(A)) (fc) dA = (Q 1 ,M>E/ ^ 2 (Ay fc (A)<A (fc) (A)<A (fc+1) (A)dA 



fe=0 



(69) = (fii^Olivn-j, 

Hence, putting all the estimates together, we conclude the proof of the lemma. 



□ 
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Lemma 9.3. Let us assume that /j, is such that L (w 2 (A) p 2 (A)) p(dX) < oo. For e W swc/i £/i<z£ 



"0 G Wg J we ftawe 

(70) 

Proof. By definition we have 



{<j>,Q 2 '^)_j = o(\\<f>f_j) 



{Q^,^)j = f3 c (C 2 ^,p)J2 w 2 (X)p 2k (X)X^^(X)dX 



k=Q 



/^/iMIx / (^ 2 (A)p 2 (A))-V(rfA) x 

\ 1/2 



1/2 



«/(A)A 2 dA ||^llo+ / w 2 (X)p 2 (X)dX\ UW, 



< cu\\ 2 j 

(71) = C\\4>\tj 

which concludes the proof of the lemma. 
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Remark 9.4. Notice that the condition on the integrability of p that appears in the statement of Lemma 
\9.3l L (ii; 2 (A)p 2 (A)) p(dX) < oo, is equivalent to assuming that p has finite moments at least up to order 

2/3 if the weights are chosen such that p(X) = y/1 + |A| 2 and w(X) = (l + |A| 2 ) and j3 < 0. 
Lemma 9.5. For cf> g W - ' 7 such that ip S Wq +1 we have 



(72) \\cr4>tj-\P S ^( x )\ 2 °{\\^tj)+\P S ^{x)\ 2 J2l w 2 {X)p 2k {X)\x^ k+1 \X) 



fe=0 



dX 



Since ip £ W Q the integral on the right hand side is well defined and bounded. 
Proof. By definition we have 



<cv,c>i mnoA 

- sup 



(73) H^'^ll = sup 

(ew J llsllj (ew 

By the definition of the operator C% we have 

./ 



ICII 



(k) 



W w 2 (X)p 2k (X)^ k \X)(CtC) {k) dX = 

J r 

= [3 s a (x)Y, / w 2 (A)p 2fc (A)^(A)(AC (1) (A))^ ; dA 

fc=0'' R + 

= (3 s a (x) J2k f w 2 {X)p 2k {X)^ k) (X)(^(X)dX + J2 [ w 2 {X)p 2k (X)X^ k \X)C {k+1) WdX 

=0 ^ R + k=0 ^ R + 

:\(3 s o-o(x)\ V / w 2 (A) j0 2fc (A) A^ fe+1 )(A) + C^(A) C (fe) (A) 



.fc=0 ' 



<2A 



< |/3 S a (x)| 



53 w 2 (A)p 2fc (A) | AV^ +1 ) (A) + CV' (fe) (A) | 

\ fc=0 



Y,w 2 (X)p 2k (X)\C^W\ 2 dX 

\ k=0 



2 s 



(74) 



< |/3 S <To(aO| 



\P S vo(x)\ 



J2 / w 2 (A)p 2fc (A)|AV>( fe+1 )(A) + CV( fc )(A)| dA, 



. £/ ™ 2 (A)p 2 *(A)|C«(A)| A d\ 

\ k=Q-' R + 



J2 w 2 {\)p 2k {\) |AV-( fc+1 )(A)| +C 2 \ipW{\)\ +2C\^ k +V(\)ipW{\) d\\\(\\ 



The first inequality follows from integration by parts on the second integral and the second and third 
inequality by Cauchy- Schwartz using Condition 17. II Therefore, by canceling the term ||CI|j an( l then taking 
the square of the resulting expression, the statement of the lemma follows. □ 

Collecting the statements of Lemmas 19.21 -Lemma 19.31 we obtain the following lemma. 
Lemma 9.6. For <j) G Wq (w, p) we have 
(75) 2 (0, g*4>}_j + HCVHlj <C(1+ \b (x)\ + \a (x)f + (Q 1)At >) U\tj ■ 

Then we are in position to prove uniqueness for the limiting stochastic evolution equation. 

Theorem 9.7. Let us assume that L (iu 2 (A)/9 2 (A)) p(dX) < oo. Then, the solution to the stochastic 
evolution equation fl6\) is unique in W (w,p). 

Proof. Equation (|6ip gives, via Lemma [ 



(76) 



11**111, e/ [2($ s ,g: Xs ^ s ) 



C I E\\<f> a \\ijds. 



£l' Xs $ t 



-j 



ds 



Therefore, by applying Gronwall inequality we get that E||$ t ||_ 7 = and the inequality 



I -(.7+2) 



C || • \\_j gives E||* t 



-(J+2) 



0, which concludes the proof of the theorem. 



< 



□ 



Appendix A. Performance of Numerical Schemes 



This appendix provides additional numerical results on the performance of the numerical schemes de- 
scribed in Sections 15.3.11 and 15.3.21 

In Figure LT1 we compare the standard error of Schemes 1 and 2 when estimating the expectation of 
f(Lj,) = {Lj, — S) + (i.e., a call option on the portfolio loss). Parameters are a time step of .005, strike 
S = .12, and K = 6 as the truncation level. Scheme 1 is performed using an optimal choice of J and M given 
by (|3"Tj) . For small /3 s , Scheme 2 outperforms Scheme 1 since the conditionally Gaussian noise dominates. 
For large j3 s , the process X dominates the dynamics and the semi- analytic calculation reduces the standard 
deviation an insufficient amount to justify the additional computational time. 

We also compare the standard error of the second-order approximation (calculated using Scheme 1) and 
direct simulation of the original finite system (|T|) in Figure [5] We use a truncation level of K = 6 and a 
time step of .005 for both asymptotic and finite systems. We report exact values in the following table for 
a total computational time of 50 seconds. For N = 500, the second-order approximation has a standard 
error one order of magnitude smaller than the finite system simulation. For N = 2, 500, the second-order 
approximation has a standard error two orders of magnitude smaller than the finite system simulation. This 
means that for the standard errors to be equal, one must invest two and four orders of magnitude more 
computational resources into the finite system simulation, respectively. 





Finite System 


Second-order Approximation 


Standard Error (N = 


= 500) 


7.8436 xl0~ 4 


8.9991 xl0~ 5 


Standard Error (N = 


= 2,500) 


1.2022 xlO" 3 


2.7449 xl0~ 5 
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Figure 7. Comparison of standard error for Scheme 1 and Scheme 2, respectively described 
in Sections r5.3.1l and l5.3.2l Parameter case is T = .5, N = 500, a — .9, a = 4, A = .2, A = .2, 
and j3 c — 1. The systematic risk X is an OU process with mean 1, reversion speed 2, and 
volatility 1. 
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